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Abstract 

The effects of the final state interaction (FSI) in semi-exclusive deep-inelastic scattering of elec- 
trons off the deuteron are analyzed paying particular attention to two extreme kinematical regions: 
i) the one where FSI effects are minimized, so that the quark distribution of bound nucleons could 
be investigated, and ii) the one where the re-interaction of the produced hadrons with the spec- 
tator nucleon is maximized, which would allow one to study the mechanism of hadronization of 
highly virtual quarks. It is demonstrated that when the recoiling spectator nucleon is detected in 
the backward hemisphere with low momentum, the effects from the FSI are negligible, whereas at 
large transverse momenta of the spectator, FSI effects are rather large. Numerical estimates show 
that the FSI corrections are sensitive to the theoretical models of the hadronization mechanism. 
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I. INTRODUCTION 



According to QCD, deep- inelastic scattering (DIS) is the process in which an incident 
electron interacts with a target quark by exchanging a gauge boson, making the quark 
highly virtual. The formation of the final, detectable hadrons, occurs after the space-time 
propagation of the created nucleon debris, with a sequence of soft and hard production 
processes. The theoretical description of these processes, which generally cannot be treated 
within perturbative QCD, requires the use of model approaches. Most of them are based 
upon the quark color string model according to which at world interval of the order of 
~ Ifm, the string which is formed by the highly virtual leading quark and the remnant 
target quarks, breaks into a hadron and another, less stretchy string. Further, at longer 
space-time intervals, this decay process iterates unless the energy of the string is too low for 
hadron production and all the final hadrons are formed. However, since the hadronization 
process can also be accompanied by gluon perturbative bremsstrahlung Q], the string model 
itself is not sufficient for a consistent treatment of hadronization. A reliable model must 
incorporate both the perturbative and the non-perturbative aspects of the hadron formation 
process. Note, that the hadronization process starts at extremely short space-time intervals, 
hence a direct experimental study of these intervals is difficult to undertake in DIS off 
a free nucleon. As a matter of fact, the final hadrons do not carry much information 
about their early stage of hadronization, and therefore only nuclear targets, which consist 
of a number of scattering centers, allows one to probe short times after the hadronization 
has started. In a nucleus, at each hadronization point one expects re-interactions of the 
produced hadrons with the nuclear constituents, so that the multiplicity of final particles 
is predicted to be reduced relative to the case of nucleon targets. Thus, by comparing the 
same DIS process off a single nucleon and off nuclear targets, information on the space- 
time structure of the hadronization process could be obtained. The theoretical model of 
hadronization developed in Ref. [3], proved to be very effective for the explanation of the 
leading hadron multiplicity ratios (nucleus to nucleon) measured at HERMES P] in semi- 
inclusive processes. It should however be pointed out that the initial stage of hadronization 
is difficult to investigate by semi-inclusive processes, where the non leading hadrons are 
strongly affected by subsequent cascade processes and therefore do not carry information 
on their formation mechanism. Recently, it has been shown that the deep inelastic semi- 
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exclusive process A{e, e', {A — 1))X, where the nucleus (A — l) is detected in coincidence with 
the scattered electron, o„g.na„y proposed to ^t.gate the .t.uctwe function of bound 

nucleons[6|, could be an effective tool to study the mechanisms of hadronization and the 
initial stage of hadronization. The simplest of the processes A{e, e', {A — 1))X, namely that 
with a deuteron target, viz 

e + D = e' + X + N (1) 

has been the object of many theoretical calculations, mainly aimed at studying the neutron 
structure function 

0-0 

, whereas its experimental investigation is planned to be performed 

at JLab llo(. 

Process has many attractive features with respect to the inclusive process '^H{e, e')X. 
As a matter of fact, it should be stressed that in spite of the fact that inclusive DIS pro- 
cesses have provided us in the past with fairly precise knowledge of parton distributions in 
hadrons, conclusive information about the origin of the EMC effect is still lacking; moreover, 
important details on the neutron structure function are unknown, which is mostly due the 
difficulties and ambiguities related to the unfolding of the neutron structure functions from 



nuclear data ll|. Semi-exclusive processes could provide, on the contrary, unique informa- 
tion on both the origin of the EMC effect, and the details of the neutron structure function; 
moreover, they can also be used as a unique tool to investigate hadronization processes. Ob- 
viously, a reliable treatment of semi-exclusive processes requires a careful treatment of the 
FSI of the nucleon debris X with the final nuclear system {A — 1). Intuitively, one expects, 
on one hand, that if the proton is detected in the backward hemisphere, FSI effects should 
not play a relevant role, so that the process could be used to investigate the bound nucleon 
structure functions; on the other hand, the effects from FSI are expected to be relevant 
in the process when the recoiling nucleon is detected in the direction perpendicular to the 
three-momentum transfer, in which case information on the hadronization mechanism could 
be obtained. In both cases, a quantitative estimate of FSI is a prerequisite for obtaining a 
reliable estimate of either nucleon structure functions or hadronization processes. In view 
of the planned experiments at JLab, a quantitative calculation of FSI effects in process (pj 
is called for. It is precisely the aim of the present paper to illustrate the results of the 
calculations of process in various kinematical conditions taking FSI into account. 

Our paper is organized as follows. In Section |n] the kinematics and the general formula 
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for the cross section are presented; the theoretical reaction mechanism for the considered 
process is illustrated in Section UTTl both in the Plane Wave Impulse Approximation (PWIA) 
()111 A|) and taking into account FSI effects (jlll B|) . In the same Section (jlll C|) the effective 
cross section describing the rescattering of the nucleon debris with the spectator nucleons is 
illustrated. Eventually, the results of the numerical calculations are presented and discussed 
in Section HVl 



II. KINEMATICS AND CROSS SECTION 



The general theoretical formalism of process (HI) can be found in several papers (see e.g. 

Q n n ^ ^ 

[7|], j8[, and p), therefore only few general aspects of the problem will be recalled here. 

In one-photon-exchange approximation, the cross section for the process can be writ- 
ten as follows 
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where a is the fine-structure constant, 
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Bjorken scaling variable, y = v/ E, the momentum of the detected recoiling nucleon (also 
called the spectator (s) nucleon), 1^^ and the electron and deuteron electromagnetic 
tensors, respectively; the former has the well known standard form, whereas the latter can 
be written as follows 



^Jl = E(P^I^m|P/)(P/I^-|Pi5)^^'^ {k + Po-k' -px-Ps)dTx, 



(4) 



where is the operator of the deuteron electromagnetic current, and Pd and P/ = Px + Ps 
denote the three- momentum of the initial deuteron and the final hadron system, respectively, 
with px being the momentum of the undetected hadronic state created by the DIS process 
on the active nucleon. In Eq. (jS} the various Wi represent the nuclear response functions 
and 
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It is well known, that within the PWIA, i.e., when FSI effects are disregarded, the 
four response functions in Eq. (jS)) can be represented only via two independent structure 
functions, viz Wl and Wt- Moreover, in the DIS kinematics, when the Callan-Gross relation 
holds, i.e. 2xFf (x) = F^{x) ( -Fj^(a;) are the nucleon DIS structure functions in the Bjorken 
limit), the semi inclusive cross section (jSJ will depend only upon one DIS structure function, 
e.g. F^{x). In the presence of FSI all four responses contribute to the cross section Q; 
however if FSI effects are not too large, nucleon momenta are sufficiently small and the 
momentum transfer large enough, one can expect that the additional two structure functions 
are small corrections, so that the semi-exclusive DIS cross section can still be described by 
„„e. effective .tn,ct„.e function Fiti<f,.,,, (BBS) (the two st.uetn.e fnnet.ons 
and Wtt can exactly be eliminated in the parallel kinematics , i.e., by choosing the spectator 
nucleon along the momentum transfer q, i.e., 6'^ = or 6*^ = vr, or integrating the cross section 
over 0, Q). 

The cross section can be also represented in terms of the light cone variable as, p\\ 
and Pt, defined as 

cts = — — —, P\\ = \Ps\ cos6'5, Pt = \Ps\ sm9s (6) 
m 

where dg is the angle between and q and the spectator four-momentum is ps = {Es,Ps) 
(note, that in the DIS kinematics the light cone z-axis is directed opposite to the vector q). 
In terms of the above variables the cross section will read as follows 
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TT^mz/ 
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III. THE SPECTATOR MECHANISM 



In what follows we will consider the reaction (pj) within the so called 
spectator mechanism, according to which the DIS process occurs on the so called active 
nucleon, e.g. nucleon "1", and the second nucleon (the spectator) recoils and is detected in 
coincidence with the scattered electron. At high values of the 3-momentum transfer q the 
produced hadron debris propagates mostly along the q direction and re-interacts with the 
spectator nucleon. The wave function of such a state can be written in the form 



*/({0)rx-,r.) = 0/3^({O)V'p.,p.(rx,r,), 
5 



(8) 



where r^; and are the coordinates of the center-of-mass of system X and the spectator 
nucleon, respectively, {^} denotes the set of the internal coordinates of system X, described 
by the internal wave function (t>i3f{{^})i 13 f denoting all quantum numbers of the final state 
X; the wave function ■?/'p^^p^(r^., r^) describes the relative motion of system X interacting 
elastically with the spectator s. The matrix elements in Eq. (0]) can easily be computed, 
provided the contribution of the two-body part of the deuteron electro-magnetic current 
can be disregarded, which means that the deuteron current can be represented as a sum 
of electromagnetic currents of individual nucleons, J^{Q'^,X) = + j/^^- Introducing 
complete sets of plane wave states |pi', P2' ) and |pi, pi ) in intermediate states, one obtains 



{Pn\j;:\Pf,Pf = p. + Ps) = 

J2 J2 (Pl3|Pl',P2')(Pl',P2'|jJ''|/5,Pl,P2)(/?,Pl,P2| /5/,p^,Ps) = 
/3,pi',P2' P1.P2 

/ (|f3^^(P)(Pb'f (Q^P■g)l/5/,P + q)^K,(p + q/2), (9) 

where = (px — Ps)/2. In Eq. Q the matrix element ( p|j^^(Q^, k-q)\Pf,p + q) describes 
the electromagnetic transition from a moving nucleon in the initial state to a final hadronic 
system X in a quantum state jSf. Obviously, the sum over all the final state /5/ of the 
square of this matrix element, complemented with the corresponding energy conservation 
5-function, generates the deep inelastic nucleon hadronic tensor of a moving nucleon. 



A. The PWIA 



In PWIA, 7* interacts with a quark of the neutron, a nucleon debris is formed and the 
proton recoils without interacting with the debris. The relative motion debris-proton is thus 
described by a plane wave 

(p + q/2) ~ (27r)35(3) (p + q/2 - a.^) = {27rf6^'^ (p - p.) (10) 

and the well known result 

'^'^ ir(x,g^p.)n^(|p|)F2^/^(Q^x,p.), (11) 



is obtained, where p = —ps is the momentum of the struck nucleon before interaction, 
K{x,Q'^,Ps) a kinematical factor (see, e.g. ref. Q), F^^'^{Q'^,x,Ps) = 2xF^^'^{Q'^,x,Ps) the 



DIS structure function of the active nucleon, and ud the momentum distribution of the hit 
nucleon, i. e. 

^i?(|p|) = 37^E y c^'r^i,A^^(r)expHpr/2) . 



Ml 



(12) 



In ffeneral, since the active nucleon inside the deuteron is off-mass shell, its deep inelastic 

n n 

electromagnetic tensor has a more complicate structure [16|, ll7[ containing off mass shell 
corrections. However, at low values of the 3-momentum of the spectator, i.e., at low 3- 
momenta of the active nucleon before interaction, these corrections are negligible and one can 
safely express the nucleon tensor only in terms of the DIS structure function F^^^{Q'^, x,Ps)- 



B. FSI 



Consider now FSI effects within the kinematics when the spectator is slow, and the 
momentum transfer large enough so that the rescattering process of the fast system X off 
the spectator s could be considered as a high-energy soft hadronic interaction. In this case 
the momentum of the detected spectator only slightly differs from the momentum p 
before rescattering, so that in the matrix element ( p\jf^^{Q^,P ■ <l)\Pf, P + q) one can take 
p ~ — Ps, obtaining in co-ordinate space 

(P^Ij^IP;) - j^(Q2,x,p,) I dV^^(r)^+/r)exp(zrq/2). (13) 

The cross section then becomes 
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K{x, Q\p,) nl^\p,, q) F2"'"(Q^ x, p.) 



dxdQ'^ dagdpj, 



N/A, 



where 
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[Vs.f\) = H J c?^'"^i,.MBW^K/r)exp(irq/2) 



(14) 



(15) 



(2VV Mo 

is the distorted momentum distribution, which coincides with the momentum distribution 
of the hit nucleon fEq.()12j)) when ~ exp(— i^jr) . 

All the effects from the FSI in Eq. fll4j) are implicitly generated by the dependence of the 
distorted momentum distribution ()15|) upon q, and by the features of the rescattering wave 
function -j/^^^ (r); the latter describes the rather complicated many-body system X+spectator 



and can only be generated by model approaches. However, in our particular case, when the 
relative momentum ~ q is rather large and the rescattering processes occur with low 
momentum transfers, the wave function ip^^ (r) can be replaced by its eikonal form describing 
the propagation of the nucleon debris formed after 7* absorption by a target quark, followed 
by hadronization processes and interactions of the newly produced pre-hadrons with the 
spectator nucleon. This series of soft interactions with the spectator can be characterized 
by an effective cross section aeff{z,Q^,x) depending upon time (or the distance z traveled 
by the system X). Thus the distorted nucleon momentum distribution (Eq. [T3j) becomes 

(16) 

where Xf is the spin function of the spectator nucleon and ^(r, q) the S'-matrix describing 
the final state interaction between the debris and the spectator, viz. 

Sir, q) = 1 - e{z) ^effi^'Q'^'^^^i^-^^) eM-bV2bl). (17) 
where the z axis is directed along q, i.e. r = z-^ + b. 

|q| 

As previously mentioned, such a treatment of FSI holds if, during the process of 
hadronization, the FSI occurs as a series of soft, elastic scattering of the produced de- 
bris with the spectator nucleon and the kinematics is chosen in such a way, that the DIS 
structure function of the nucleon does not change too rapidly. In this case the factorization 
assumption holds and the rescattering process is described as the debris-nucleon rescattering 
via an effective cross section. Hence, the usual eikonal approximation can be applied. Note, 
that the above expressions can also be used to calculate the effects of FSI in quasi-elastic 
scattering (QES) by simply replacing the debris-nucleon cross section with the z-independent 
nucleon-nucleon cross section (see, e.g., [l3,I3])- 



n 



FSI I 

D \ 



i i f 

P.,q) = 37^^I] J c^r^i,AiB(r)'5(r,q)x/ exp(-ip,r) 



C. The effective debris-nucleon cross section 

Recently in Ref. the final state interaction in DIS off nuclei due to the propagation of the 
struck nucleon debris and its hadronization in the nuclear environment, has been considered 
and applied to the process of the type A{e, e'{A — 1))X with A> 4. This process precisely 
coincides with the one investigated in the present paper, i.e the process ^H{e, e'p)X. The 
basic ingredient governing the FSI, i.e. the time- and Q^-dependent effective cross section 
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(Te// describing the interaction of the debris with a nucleon of the spectator system (A — l), 
has been obtained in Ref. on the basis of a model which takes into account both the 
production of hadrons due to the breaking of the color string, which is formed after a quark 
is knocked out off a bound nucleon, as well as the production of hadrons originating from 
gluon radiation. The general expression has the following form 



NN 
tot 



ttN 
tot 



nM(t) + ncit) 



(18) 



where cr^/^ and are the total cross sections of nucleon-nucleon and meson-nucleon 
interaction, and UMif) and ncit) are the effective numbers of created mesons and radiated 
gluons, respectively, and are explicitly given in Ref. 0. There the color-dipole picture was 
employed by replacing each radiated gluon by a color-octet qq pair. 

The cross section (|18p exhibits a rather complex Q"^- and x-dependence, which, however, 
asymptotically tends to a simple logarithmical behavior. This is illustrated in Fig. ^ where 
the Q^-dependence of the effective cross section, calculated with the set of parameters given 
in Ref. 5|, is exhibited. 



IV. RESULTS AND DISCUSSION 

Eq. ()14p is the basic equation that one has to evaluate in order to provide a quantitative 
significant description of the quasi-exclusive process A{e,e'{A— 1))X; the latter, as already 
mentioned, has been originally suggested l(| as a powerful tool to investigate the reaction 
mechanism of DIS off nuclei and the bound nucleon structure function; it was however soon 
realized P| that it could also be an extremely powerful tool to investigate the hadronization 
mechanism. In the present paper we will focus on both aspects of the problem, namely we 
will address the problem of finding two different kinematics, a one in which FSI corrections 
are negligibly small, so that a direct study of the DIS nucleon (neutron) structure functions 
becomes possible, and another one, where FSI are maximized, thus obtaining information 
on the hadronization mechanism. As it can be seen from Eq. ()14|1 . all of the FSI effects 
are contained in the distorted momentum distribution n^^\ps,q) (Eq. (fT^ ). so that the 
investigation of its deviation from the nucleon momentum distributions (?2/)(|p|))) could 
provide clear signature on FSI effects. This is illustrated in Fig. |21 which shows the ratio 
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^rf^/^D, (with n^'^^ and no given by Eqs. (fT^ and (fT^ . respectively) calculated using two 
different models for the effective cross section o"e//, representing its upper and lower limits, 
viz the time- and Q^- dependent cross section of Ref. 5| (solid lines), and a constant cross 
section a^ff = 20 mb (dashed lines) considered in Ref. 8|. In our numerical calculations the 
parameters entering Eq. (|17p. viz the slope bo and the ratio a of the real to the imaginary 
part of the forward amplitude, have been taken from A^A^ scattering data at high energies. It 
should be pointed out, in this respect, that our results, in the range of considered momenta, 
are not very sensitive to the value of a; whereas the values of the latter is known in the 
case of nucleon-nucleon scattering, the value for debris-nucleon scattering should rely on 
some theoretical models. This point is under investigations and the results will be presented 
elsewhere. It can be seen from Fig. |21 that the predictions given by the two different models 
of the effective cross section ere rather different, particularly when the recoiling proton is 
emitted perpendicularly to q, i.e. at 6^ ~ 90°, and with large values of the momentum {ps ~ 
0.2 GeV/c). Thus, by investigating this kinematical region one could, in principle, obtain 
unique information about the magnitude of (Je// and, consequently, about the hadronization 
mechanism. 

In order to illustrate the different role played by the FSI in Deep Inelastic and Quasi- 
Elastic scattering, we show in Fig. Elthe ratio n^^^ /n£, for the QE process ^i/(e, e'n)p, which 
corresponds to the calculation of Un^^ using the effective nucleon-nucleon cross section, i.e. 

n 

(Te// = 44 mb. It can be seen, in agreement with the results for heavier nuclei p], that the 
survival probability of {A — 1) in DIS is less than in QE scattering. 

The results exhibited in Figs. |21 and IHl also demonstrate that FSI effects are essentially 
reduced in parallel kinematics {6 = 0°, 180°) and also at small values of Ps- This region 
is attractive from the experimental point of view since in such a kinematics the structure 
functions Wlt and Wtt in Eq. Q vanish exactly and the effective structure function F2 in 
Eq. (|TT| exactly coincides with the DIS structure function of a bound nucleon. Thus, in this 
kinematical region a direct experimental study of the off-mass shell effects in DIS becomes 
possible. The results exhibited in Figs. El and El have been obtained by u sing the deuteron 
wave function corresponding to the Reid Soft Core (RSC) potential 1^. At small values 
of IpI different potential models provide basically the same deuteron wave function, which 
is not true, however, at moderate and large values of |p|. Therefore, we have investigated 
the dependence of the ratio n^^^ /no upon the two body potential used to generate the 
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deuteron wave function. The results are presented in Figs. |3] where the predictions of the 
Reid and Bonn potentials are compared in correspondence of two different values of the 
transverse momentum pt (cf.Eq. IB)). It can be seen that, independently of the emission angle 
of the spectator nucleon, the two potentials provide essentially the same results at relatively 
small values of the spectator momenta \ps\, but very different results at large values of \ps\ 
{\Ps\ > 0.25GeV/c). 

The dependence of the ratio n^^^/riD upon the electron kinematics (Q^ and x) is en- 
tirely contained in the and x-dependence of the effective cross section cTg// given by Eq. 
(fll^jl . Since this dependence is weak (see Fig. Hj), FSI effects in the quasi-exclusive process 
'^H{e,e'p)X are practically and x-independent (see, e.g. Figs. |S]and|ni)- 

To sum up, from the analysis we have exhibited, it can be concluded that FSI effects in 
the semi-exclusive process are negligible in the backward kinematics with slow momenta 
of the detected nucleon, which would allow one to investigate the nucleon structure function 
of bound nucleons, in particular the neutron one; if, on the contrary, the spectator nucleon 
is emitted perpendicularly to the momentum transfer, FSI effects are enhanced and different 
models for the hadronization process could be investigated. 

This work was partially supported by the Ministero dell'Istruzione, Universita e Ricerca 
(MIUR), through the funds COFINOl. L.P.K. is indebted to the University of Perugia and 
INFN, Sezione di Perugia, for warm hospitality and financial support. 
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FIG. 1: The debris- nucleon effective cross section (Eq. I18() plotted vs the distance z for a fixed 
value of the Bjorken scaling variable x and various values of the four-momentum transfer (after 
Ref. 0) 
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FIG. 2: Left panel: the Deep Inelastic Scattering ratio rij^^^ /ud with Uj^^^ and ud given, respec- 
tively, by Eqs. (fT^ and (fT^ . calculated vs. the momentum ps = \ps\ of the spectator nucleon 
emitted at different angles 9s- The full lines correspond to the Q^- and z-dependent debris- nucleon 
effective cross section (Jeff shown in Fig. whereas the dashed lines correspond to a constant 
cross section a^fj = 20 mb. Right panel: the same as in the left panel, vs the spectator emission 
angle 9s for different values of the spectator momentum. Calculations have been performed at 
= 5 {GeV/cf and x = 0.2. 
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FIG. 3: Left panel: the Quasi-Elastic Scattering ratio n^'^^/n^) with Uj^^^ and ud given, respec- 
tively, by Eqs. (fTO)) and (fT^ . calculated with the nucleon-nucleon (Jeff, i-e. a^ff = 44 mb, vs. 
the momentum ps = \ps\ of the spectator nucleon emitted at different angles 6s ■ Right panel: the 
same as in the left panel, vs the spectator emission angle 6s for different values of the spectator 
momentum. Calculations have been performed at = 5 (GeV/c)'^ and x = 0.2 . 
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FIG. 4: The Deep Inelastic Scattering ratio juY) vs the light cone variable Os (cf. EqlHl), 
calculated at = 5 GeV'^ /(? and x = 0.2 using two different deuteron wave functions: the solid 
lines correspond to the RSC potential and the dashed lines to the Bonn potential. Calculations 
have performed in correspondence of two values of the transverse momentum px of the recoiling 
nucleon (cf. Eq. namely = GeV/c (open dots) and px = 0.2 GeV/c (full dots). The results 
correspond to the debris-nucleon cross section shown in Fig.l. 
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FIG. 5: The Deep Inelastic Scattering ratio juY) vs the light cone variable Os (cf. EqlHl), 
calculated at different values of and fixed values of x = 0.2 and pT = GeV/c. The deuteron 
wave function corresponds to the RSC interaction. 
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